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ABSTRACT 

Gravitational waves (gw) propagating through a magnetised plasma excite low- 
frequency magnetohydrodynamic (mhd) waves. In this paper we investigate whether 
these waves can produce observable radio emission at higher frequencies by scattering 
on an an-isotropic intrinsically relativistic distribution of electrons and positrons in 
the force-free wind surrounding a double neutron star binary merger. The relativistic 
particle distribution is assumed to be strictly along the magnetic field lines, while 
the magneto-plasma streams out at a relativistic speed from the neutron stars. In 
the case of Compton scattering of an incident MHD wave transverse to the magnetic 
field, we find that the probability of scattering to both a transverse x-mode and a 
quasi-longitudinal Langmuir-o mode is suppressed when the frequency of the incident 
wave is below the local relativistic gyro-frequency, i.e. when the magnetic field is very 
strong. 

Key words: (magnetohydrodynamics) MHD - gravitational waves - magnetic fields 
- plasmas - radiation mechanisms: general - scattering 



1 INTRODUCTION 

This decade is expected to witness the historical first direct 
detection of gravitational waves (gw) with detectors such 
as the Advanced Laser Interferometer Gravitational Wave 
Observatory (ligo 1 ) and the Laser Interferometer Space 
Antenna (LISA 2 ). The former is designed to detect gravita- 
tional waves from compact sources such as merging neutron 
star binaries in the kHz regime, whereas the latter will de- 
tect the lower frequency (< 1 Hz) gravitati onal waves from 
thous ands of compact white dwarf binaries l)Nelemans et alj 
l200ll) as well as the spiral-in signal of su per-massive bl ack 
holes, throughou t the visible universe fsee lSchutzl l|l99fll ) or 
iKokkotasl ll2004) for a review of likely sources for LIGO and 
LISA). 

Gravitational waves are emitted by highly energetic 
events but occur at relatively large distances and the sig- 
nal that reaches Earth is exceedingly weak. To extract the 
GW from the noisy signal, some theoretical knowledge of the 
expected waveforms is essential. To confirm the detection of 
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a GW burst, any additional electromagnetic signature of the 
event would be extremely useful. 

It so happens that many of the possible GW sources are 
embedded in a strong magnetic field. Examples are rapidly 
spinning neutron stars which have a small oblateness and 
precess, accrete, or are unstable to the excitation of the r- 
mode, supernova core collapses and bounces, newly born 
'boiling' and oscillating neutron stars, magnetars with crust 
fracturing (~ 20 Hz), and coalescing compact binaries in 
which at least one component is a magnetic neutron star. 

In the last case, numerical models pred ict that maxi- 
mum GW luminosities of the order of 10 48 W jlbrahim et alJ 
2003) are released into a wound-up magnetic field of up 
to 10 8 - 10 12 T. We have investigated i n previous papers 
jMoorteat fc Kuiiperd feOOl 12004 120051) referred to from 
hereon as Paper I — III) whether these extreme space-time 
distortions perturb the ambient magnetic field sufficiently 
to produce an observable electromagnetic counterpart of the 
GW burst. We have found that a GW can couple linearly at 
the same frequency to all three low-frequency plasma eigen- 
modes. The shear Alfven mode and the transverse slow and 
fast magnetosonic waves (msw) are excited depending on 
the polarisation of the gravitational waves with respect to 
the ambient magnetic field orientation, but only the fast 
magneto-acoustic mode can interact coherently with the GW 
a over longer time (distance) and is therefore the most in- 
teresting. 

The MSW can not escape directly as observable radia- 
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tion. In Paper III we pointed out that the MSW encounters 
relativistic electrons and positrons with Lorentz factors dif- 
ferent from the bulk Lorentz factor of the particles carrying 
the MSW, for instance due to a tail in the distribution func- 
tion or simply due to inhomogeneities in the plasma such 
as beams. For typical Lorentz factors of a few hundreds, 
the MSW would than be inverse Compton (ic) scattered to 
low-frequency radio emission in the LOFAR (the LOw Fre- 
quency Radio Array 3 ) band. In a first estimate, we approxi- 
mated the scattering cross-section with the Thomson cross- 
section and predicted a flux that would be easily detectable 
by LOFAR for a double neutron star merger at a distance 
of ^ 1 Gpc. Clearly, the indirect detection of gravitational 
waves with a radio array would be of extreme interest. 

In the present paper we study in more detail the scat- 
tering process of a fast MSW into escaping electromag- 
netic waves in a strongly magnetised an-isotropic relativistic 
plasma. 

The paper is set up as follows. In Section [5] we reca- 
pitulate the intuitive picture presented in Paper II which 
underlies the excitation of mhd waves by a gravitational 
wave propagating through a magnetised plasma. To in- 
vestigate the subsequent scattering of these mhd waves 
we need a covariant relativistic theory for p lasma dynam- 
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emission processes in pulsar magnetospheres. We apply this 
theory to the force-free wind well outside the effective light 
cylinder of the merging binary. In particular, we consider a 
plasma which -in the frame comoving with the wind- has 
a momentum dispersion along the magnetic field only, and 
construct the covariant equations for the GW-excited cur- 
rent, the linear and non-linear response to that current and 
the scattering probability in Sections in Section |Bj we 
apply these general results to the geometry of an incident 
transverse magnetosonic wave that scatters on electrons and 
positrons with a 1-dimensional relativistic dispersion along 
the magnetic field and calculate the probability that either 
transverse or longitudinal waves are emitted in the direc- 
tion of the observer. A numerical analysis is presented in 
Section [7] where we calculate the energy radiated into scat- 
tered radio waves. The conclusions follow in Section E] 

Standard geometrised units are used throughout this 
discussion (G = c = 1), Greek indices are used for time- 
space components (fi = . . . 3) and Latin indices for spatial 
only components (i = . . . 3). 




Figure 1. Deformation of the world-sheet of a ring of test parti- 
cles by the tidal force of a passing + polarised GW (with k gw _LB°). 
The magnetic field is stretched and compressed and excites a 
magneto-compressional wave. 



In this section we discuss what happens when the two neu- 
tron stars at the origin of the system merge and release a 
large amount of GW into the wind. 

2.1 Gravitational wave tidal field 

The driving force exerted by a GW on test particles is de- 
scribed by Einstein's field equations or alternatively by the 
general relativistic equations for geodes ic deviations de scrib- 
ing tidal accelerations (see for instance IStewartl jl990t0 . In- 
tegrating those differential equations twice results in the 
well-known spatial deviations of test masses in an interfer- 
ometer detector such as LIGO: 



5x = -(h+x + h x y ) 



x 



h+yo) , (1) 



where xo and yo are the initial separations of test-masses 
along the x- and y-axis in an L-shaped interferometer, 8x 
and 8y are the excursions causes by the passing GW and h+ 
and hx refer to the two possible polarisations of the GW. 

2.2 Interaction with magnetic field 



2 GW => MAGNETO-ACOUSTIC WAVE 

In a pulsar environment plasma flows out along the open 
field lines and develops into a force-free wind outside the 
light-cylinder in which the toroidal component of the mag- 
netic field quickly dominates the poloidal one (see Section 
and Fig. . The magnetic field is therefore essentially per- 
pendicular to the radial propagation of the wind as in Fig.0 



LOFAR website: |http:/ /www . lof ar . org/ 



In the ideal mhd approximation the magnetic field lines are 
frozen into the plasma, and the magnetic field will exhibit 
the same excursions as the test particles. The action of a GW 
propagating in the z direction is only in the x — y plane, and 
when we choose the orientation of the ambient magnetic field 
in the x — z plane, a GW which is + polarised with respect 
to these axes (so h x — 0) excites a growing magnetic field 
perturbation with: 



SB X oc \h+B° x 



.W„ x \h+Bl = , (2) 



The resulting interaction is illustrated in Fig. Q where we 
show the world-sheet of a ring of test particles in the tidal 
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Figure 2. The orientation of the MSW components and the per- 
turbed magnetic vector field (inset). 



field of a passing GW, which we choose to propagate in the z- 
direction (and k gw _LB°). Since the GW amplitude depends 
only on z — t, we can interpret the vertical axis either as 
the propagation direction 2 at a certain time to or as the 
temporal evolution for a ring of particles around the spatial 
origin z = zq as in a Minkowski diagram. In the absence 
of the GW the world-sheet would be a straight cylinder, but 
due to the GW the circular circumferences are periodically 
stretched and compressed into ellipses with axes along the 
x and y axis for this particular + GW polarisation. The cor- 
responding behaviour of a magnetic field that was spatially 
uniform in the absence of the GW, is a periodic amplifica- 
tion of the magnetic field when B® is amplified by the GW 
(Bx > 0) and dilution when B® is suppressed (B* < 0). The 
result is a magneto-compressional MSW mode. 

2.3 Magnetosonic waves in a Poynting flux wind 

When the plasma is tenuous and strongly magnetised such 
that the magnetic pressure greatly exceeds the gas pres- 
sure, the sound velocity is much smaller that the relativistic 
Alfven velocity which approaches the speed of light ua — c. 
In this Poynting flux dominated limit we can neglect the gas 
pressure and the MSW and the GW obey almost the same 
dispersion relation tj gw = w msw = fc gw c = i msw MA and can 
interact coherently over the longest interaction scales. 

The orientation of the different components of the fast 
MSW are shown in Fig. The polarisation is transverse, 
such that the electromagnetic field components and the wave 
vector form a mutually orthogonal triad: E 1 _LB 1 _Lk. Also 
shown is a schematic illustration of the perturbed magnetic 
field for on oblique magneto-acoustic wave propagating in 
the 2-direction. The perturbations are exaggerated to em- 
phasise the effect. Since B z is constant while B x oscillates, 
the total magnetic field has an overall wavy pattern. 

As was mentioned in the Introduction, the magne- 
tosonic waves excited by a GW have frequencies in the kHz 
regime and cannot escape the plasma directly. However, 
upon inverse Compton scattering of the MSW on relativis- 




Figure 3. Merging neutron star binary 



tic particles with a typical Lorentz factor 7 of a few hun- 
dreds, the photon frequency can be boosted by up to a fac- 
tor 27 2 larger than that of the incident MSW wave and into 
the LOFAR band (30-240 MHz). Close to the source, such 
frequencies are still well below the local plasma frequency 
as determined by the pair plasma expelled by the source 
which effectively behaves as a single millisecond pulsar with 
a strong magnetic field. The region where scattering into es- 
caping radiation is most likely to occur is far away from the 
pulsar where the plasma frequency has dropped sufficiently 
(see also Section 0. 

In the following sections we will first summarise the co- 
variant description of some fundamental ingredients of rel- 
ativistic plasma dynamics and then derive the probabilities 
of scattering a MSW into the higher frequency normal modes 
of an intrinsically relativistic plasma. Note that while in the 
context of this section it was most illustrative to choose the 
GW wave vector along the z axis and the magnetic field at 
some angle to it, in the remainder of this paper it will be 
more convenient to do the opposite and choose B° along the 
z axis and the wave vectors in arbitrary directions. 



3 EMISSION 

3.1 Wave equation and photon propagator 

We will study the emission due to the external current gen- 
erated by a passing mhd wave. The inhomogeneous wave 
equation for the scattered field of arbitrary wave mode P 
with vector potential ^4p is given by: 

A^(fc)A Pl ,(fc) = -^o-CtW , (3) 

where is the wave operator for the homogeneous equa- 
tion which determines the polarisation properties of wave 
mode P and includes the induced self-consistent field. In 
the weak turbulence regime one can expand the current ex- 
cited by a field J 4 M (fc) in terms of a power series in this field. 
The first two terms in this expansion define the linear and 
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quadratic non-linear response: 

x (27r)V ra (-*,fc,fa)A"(Jfei)A a (fe) , (4) 

where stands for the 4- vector (u>i, kj). A plasma theory is 
different from a vacuum theory in that the linear response 
of the plasma is included in the wave operator: 

A M "(fc) = fcV - k"k v + noa^ik) , (5) 

which follows from Maxwell's equations ik v F^ v = /io J M with 
F^ LV = ik^A" - The homogenous wave equation (JSJ 

with the external current set to zero), and © for the wave 
operator then determine the wave eigenmodes of the plasma. 

In Section ^] we will calculate the external current ex- 
cited by a (GW induced) magnetosonic wave propagating 
through a magnetised relativistic plasma and in Section 
we will derive the response of the plasma to this current. 
To find the associated scattered wave field, has to be 
inverted: 



A£(fc) = -D%{k)J^ t {k) 



(6) 



Z) M " is the Green's function for the inhomogeneous wave 
equation or the photon propagator. We are interested in the 
radiation field excited by the external current so we only 
consider the irreversible anti-Hermitian part of the photon 
propagator which is obtained by imposing the causal condi- 
tion on it and using the Plemelj formula. 

The covariant and gauge independent expressions for the 
photon propagator and its anti-Hermitian part are 4 : 

G a G'p X^(k) 



D""{k) 



I Hi 



D 



(*) 



(Gk)(G'k) 

Ga.G 



X(k) 



-X 



fj,auf3 



(k)5[X(k)] 



(7) 



(8) 



(Gk)(G'k) 

where X^ avli and X(k) are the second-order matrix of cofac- 
tors and the determinant of A M ", respectively. G M determines 
the gauge and can be chosen, for instance, as the Lorentz 
gauge Gfj, = Coulomb gauge G M = [0, k] or temporal 
gauge G M = [1,0]. 

We will specify to the temporal gauge as it is the most 
convenient and the only gauge that allows a unique nor- 
malisation of the wave polarisation 4-vectors: one is free to 
specify e^ M (k)e* M ^(k) — —1 with e^ M — [0,ejj]. Furthermore, 
in the temporal gauge the first-order matrix of cofactors and 



the determinant of the 3-tensor A* 3 - are given by A' 



A u 



0) 



and A (t) (fc) = Lu 2 X{k), respectively, with A'jA'j- = A (t) ^. 
The polarisation vectors can be derived from the spatial ma- 
trix of cofactors by 



e M i{k)eMj{k) = — 



(9) 



The photon propagator only has contributions from the 
poles at the zeros of A(&m) corresponding to the dispersion 
relation k = fcjwr for a specific wave mode. To make this more 
explicit we use the identity: 

5[w-WM(k)]+<5[w+WM(— k)] 



5[X(k)] = Y, ■ 



M 



\ax(k)/du\ 



(10) 



which can be written in terms of the ratio of electric to total 
energy 



R M (k) 



W[f\k) 
W M (k) 



a\(k) 



X\(k) 



T L 



(11) 

Here the wave energy W M (k) = P%f(k) = T$(fc) is the tem- 
poral component of the 4- momentum P^ik) and the time- 
time component of the energy-momentum tensor T^(k). 

Combining l|%}l- l)llft . the radiation field in mode M is 
found by retaining only the resonant anti-Hermitian part of 
the photon propagator, given by: 



+ e*M(k)eM(k)5[iV + ui M (-k)]} 



(12) 



while D Atlv (k) = J2m D tr ( k ) and since CH is derived us- 
ing the temporal gauge, only the spatial components (D^iYj 
are non-zero. 



3.2 Emission probability and 4- momentum 
radiated 

The 4-momentum radiated in an arbitrary wave mode due 
to an arbitrary source, described by an extraneous current, 
is identified by the work done by this current on the elec- 
tromagnetic field. Written covariantly and averaged over a 
truncation or normalisation time T (which we can let go to 
infinity) , using the power theorem and the charge continuity 
equation k" J„ = 0, one has: 



i Jd 4 x JZ*(x)FS(x) = 



(13) 



i J -0s n-ik»j: Kt (k)A„(k)] = J 



tTk 



r(*) 



Q%j(k) = dP^j(k)/dt is the rate at which 4- momentum is 
radiated into mode M. Using @ with (II l2fi and performing 
the integral over the delta functions we find: 



HoRM(k) 
Tujm 

= k^jWM( k ) 



kl, \e* Mv {k)J^(k M )Y 



(14) 



Equation 114|l defines the probability WM(k) per unit time 
and unit volume of k-space of spontaneous emission of a 
wave quantum. 

The rate of emission and the ratio of electric to 
total energy Rm can also be written more conveniently in 
terms of the eigenvalues of the response tensor 



such that 
1 



om (fc) = eM M (fc)eM 1 /(fc)a'" / (fc) 



— -^-u M {k) 



R M (k) 

M . . 2i^ R M {k)W M {k) ui A fi 



(15) 

(16a) 
(16b) 



4 In J^J and throughout this paper, inner products are denoted 
as Gk = g^G^k". 



where the purely imaginary anti-Hermitian part of l|15|l . 
a^f(fcM), measures the effect of dissipation in mode M. 
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4 GW INDUCED CURRENT 
4.1 Zeroth order current 

The zeroth-order motion of electrons and positrons in a uni- 
form and static magnetic field is a spiralling orbit around the 
magnetic field given by the covariant equation of motion: 

du"(r) 



Iu{t) 



(17) 



dr m 

where q is the charge of the particle, m the mass and r the 
proper time. One can define a magnetostatic field as having 
F^Fpv > in any frame in which case it is possible to 
choose a particular frame where E = and B/0. We will 
orient the axes such that the magnetic field lies along the 
z-axis and define B = \ \J Fg" Fo^. If we specify the initial 
condition for the 4-velocity as Uq = 7(1, v±, 0, Vv, ) at r = 
we can integrate Ea. 1 171 twice to find the orbit 5 : 

U<°>"(T) = i^(T)u 0v 



(18a) 

x£ + t^(T)u „ , (18b) 
where t^" (r) and its Fourier transform r^ v (iv) are given by: 



t M "(r) 



» = 



1 

I 1 




\ 










V 









— sin Q c t 

-J) COS Q C T 









r\ cos Sl c r 
— sin Q c r 




\ 











(19) 



-1 / 



7, vj_, Vu and the components of the 4- momentum p M = 
m/, p± — "/mv± and pu = jmvu are all constants of 
the motion. Also, we have defined r\ = qj\q\ and the (non- 
relativistic) gyro- frequency Q, c — \qB/m\. 

Note that in the limit £l c — > 0, r A " y (tj) reduces to the 
Minkowskian metric tensor g^ iv = rf" v and l|18|l gives recti- 
linear motion. 

The zeroth-order single particle current in the absence 
of any perturbations is now given in Fourier space by: 



J? p (k) = qJ dr «"(T)e ifcx(r) 



(20) 



4.2 First-order current 



We now study a particle that is orbiting in a static mag- 
netic field when an mhd wave with amplitude A^(k) passes 
and perturbs the motion, by including the covariant Lorentz 
acceleration in 11711 : 



(21a) 



iq f d 4 k! _ lfclX (r) 



m J (2tt) 4 

x kiu{T)G" v \ki,u{T))Au{ki) . (21b) 
where the useful projection tensor G M " is defined by: 



ku 



(22) 



5 Note that Uq denotes an initial condition, whereas means 
the unperturbed 4-velocity and similarly for X. 



The orbits can again be obtained by integrating the first- 
order equation of motion 12 U : 



Jo 



(23a) 



I (1) "(r) = [ dr" f dr' i^(r"-r')Sl 0) (r') , (23b) 
Jo Jo 

where Si°\r) is given by inserting the zeroth-order orbit 
Eq.EBlin Ea.l2Tbl 

The first-order expansion for the current density is: 

J^(k)=qjdrj^(r)e^ m ^, (24) 



with: 



j (1)M (r) = u (1) "(r)+ikX W (r)u (0) '"{r) . (25) 



Inserting I25|l in JUJ, partially integrating with u' 1 - )M 
dX m /dr and using 1231 . results in: 



J s ( p )M (fc) = iq j dr ku (0) {r)G c " 1 {k,u (0) (r))X l 



(1) ikX 

y 'e 



(o) 



-i- Idr 



dr'' 



dr' 



d 4 ki 



(26) 



(27T)* 

X A kX<, ^-^ x0 ^M^{r" -r^A^k,) 

X kv (r)G c " i {k,u (r))k 1 u a {r')G l3v {k 1 ,v {r')) . 

Eq. I|26|l is the current associated with a single particle. To 
calculate the response tensor for the magnetised plasma one 
can use either the forward-scattering method were all the 
perturbations are included in the orbits of the particles and 
an average is made over the initial conditions, or the Vlasov 
method in which the perturbations are included in the dis- 
tribution function. In an unmagnetised plasma, the simplest 
method is to calculate the response tensor for a cold plasma 
fluid and generalise this result to an arbitrary distribution 
function by using Lorentz transformations. This method is 
generally not suitable for a magnetised plasma since there is 
no inertial frame associated with the gyrating particles. An 
exception is the case where the particles move strictly along 
the magnetic field. Since we want to study the effect of a 
mhd wave propagating through the force-free plasma wind 
relatively close to the source where the magnetic field is very 
strong and any transverse momentum is quickly synchrotron 
radiated, a strictly parallel particle distribution should be a 
reasonable assumption. 

One advantage of the small gyro-radius approximation 
is that while l|2fifl generally has to be expanded in Bessel 
functions (J a (k±R) with R the gyro-radius) to carry out 
the integrals, in the strictly parallel case the arguments of 
all the Bessel functions vanish for v±_ oc R = which leads 
to J s (0) = for all integers s ^ and Jo(0) = 1, and 



Jsp 11 (k) 



d 4 kt 



-*(fc-fci)zo 27r(5 [( fcl| -k n )u) 



m J (2tt) 4 

x G atl {k,u)T a ^k\\u)G pu {kx,u)Amu{k{) , (27) 



where A^ M denotes the incident fast magneto-acoustic mhd 
wave and since u is strictly along the magnetic field we im- 
plicitly have ku — (ku)^ — "/(uj — fc|| u ||), v ± = 0. 
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5 LINEAR AND NON-LINEAR RESPONSE 

To find the linear response tensor for the magnetised plasma, 
we use the forward-scattering method and average the first- 
order single-particle current (1261 over the distribution of 
particles F(xo,po). The distribution function is defined co- 
variantly in terms of the number cLV of world lines -one 
per particle- in an 8-dimensional phase space threading a 
7-dimensional surface element d 4 xod 4 po/dr as determined 
by: dJV dr = F(xo,po)d 4 xod 4 po- With this identification one 
can replace the integral over dr in (12611 with an integral over 
dJV dr -that is over d 4 xod 4 po weighed by F(xo,Po)~ to av- 
erage over the ensemble of particles. 

When the particle distribution is uniform in space and 
time, F(po) does not depend on the initial conditions xo 
and the current fl2( ill only depends on xo through the un- 
perturbed orbit (I18bl . Explicitly, this results in a factor 
exp[i(fc — ki)xo] as in l|27^ and the integral over xq gives 
a delta function (27r) 4 <5 4 (fc — fei) over which the ki integral 
is performed. 

An alternative and easier procedure is to use J23 which 
is valid in the small gyro-radius approximation, realise that 
2ir8 [(fc|| — fei||)w] = J dr exp [i(k\\ — ki y )ut] and replace 
this integral with the integral over phase space. The result 
is: 



'(k) = -— I ' d 4 p F{p)G ail {k,u)T afj (ku)G r)l '{k,u) 
m J 



q_ 

m 

= f^9(p^(k,u) 



7|| 



(28) 



where the last equality defines a A "'(fc, u) and the distribution 
F(p) is chosen parallel to the magnetic field, so the spatial 
part can be written as: /(p) = S 2 (p±)g(p^). More specifi- 
cally, a relativistic thermal distribution strictly parallel to 
the magnetic field is given by the one-dimensional Juttner 
distribution: 



F(p) = 2nm8(p 2 — m 2 )S 2 (p±)g(pu) 

m& T 

9(7) 



(29a) 
(29b) 



2nmK 1 {m/T) ' 

where K v are the modified Bessel or Macdonald functions, 
T is the temperature in energy units, u is the 4-velocity of 
the rest frame in which u — [1, 0] such that pu = my. 

The Juttner distribution 1291 is normalised by the (non 
invariant) number density n evaluated in the rest frame: 



d 4 p u°F(p) 



rather than the invariant proper number density 

3(7) 



pr = / d 4 p F(p) = 



dp,. 



7 



(30a) 



which does not correspond to the actual number density in 
any frame (when the plasma is not cold). 

In this special case of a strictly parallel distribution l|28|l 
can also be obtained from a fluid description as in the un- 
magnetised case, as was mentioned in the previous section. 
The covariant linear response tensor for a cold magnetised 
plasma-fluid with rest frame 4-velocity u and proper number 
density n pr is: 



i(*) = 



q 2 n pl 



G ail {k,u)T a 0{ku)G 0v {k,u) , (31) 



and the response tensor 1281 for an arbitrary but strictly par- 
allel distribution is found by interpreting the proper 
number density d 4 pF(p), which implies implicitly that for 
each particle with 4-velocity u the response is calculated in 
its rest frame using 1311 and subsequently Lorentz boosted 
along the magnetic field to the frame of interest and then 
the contributions of all particles are summed by integrating 
over d 4 p. 

The cold plasma response (1311 can be recovered from 
l|28fl in the limit F(p) = nS 4 (pu — mu) where n = n pT . 



5.1 Linear response of cold magnetised plasma 



The explicit spatial components of 1311 will be useful later 
on. In terms of the total plasma frequency lo 2 = lu 2 + + 
u>p- , the positron plasma frequency LUp + = e 2 n+ / (eom) (and 
similarly for the electron plasma frequency) and the average 
charge number £ = (n+ — n_ )/ (n+ + n_ ) measuring the 
excess of posit rons (with density n+) to electrons, one has 
jMelroselll986T) : 



p 
Mo 











(32) 



We are interested in a pure pair plasma with £ = in which 
case the linear cold plasma response tensor l|32|l becomes 
diagonal and the contributions of the electrons and positrons 
add in the (linear) Thomson scattering process, in contrast 
to for instance an electron gas which has f = — 1. 

More generally, the plasma might have some charge ex- 
cess, for instance related to the Goldreich- Julian charge den- 
sity ncj such that £ ~ naj/n to t = 1/M. M is the multi- 
plicity of secondary particles created by a primary particle 
that emits curvature radiation while it flows out along an 
open field line anchored to the polar cap. Th e curvature pho- 
tons then produce a cascade of pair creation ijGurevich et alJ 
1993). The multiplicity is largely unknown but varies from 
10 2 -10 6 . The gyrotropic terms only contribute quadrati- 
cally in the dispersion relation for the plasma, i.e. of or- 
der 1/M 2 , wh i ch is t he motivation to neglect those terms in 
iGedalin et alJ al998l) . Section V A. For very low- frequency 
waves and low multiplicity the gyrotropic terms become im- 
portant when 1/M ~ £ ~ lj/Q. c as can be seen by com- 
paring the components of 1321 . Although this might lead to 
interesting behaviour, we will assume for simplicity that the 
plasma is sufficiently neutral that £ can be neglected. 



5.2 Linear response of relativistic magnetised pair 
plasma 

In evaluating 128^ for more complicated response tensors it 
is convenient to define the average, (V), of any variable V 



by: 



(33) 



The frame in which the bulk momentum vanishes, {"/v) — 0, 
is defined as the rest frame and the non-gyrotropic compo- 
nents can be written in terms of the Alfven speed, va, a 
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mean square speed, Av 2 , 

WA = 



A, 2 =<^> 



(34) 



/<7>w p ' (7> 
and the relativistic plasma dispersion function for the phase 
velocity z = w/k\\ parallel to the magnetic field: 

1 \ r dp dg(p) 



W{z) = 



(35) 



, 7 3 (z — v) 2 I 7_oo v — z dp 

Since we are interested in scattering of kHz waves to 
~ fOO MHz waves that both have ku Q c , we can expand 
in ku/Clc and only keep terms up to second order. Under 
these approximations, the linear response tensor is given in 
terms of the refractive index n = k/uj (and n± = nsin(9, 
nil = n cos 8) by: 

■2 ( 1+njjAw 2 -n±n\\Av 2 \ 

l + njjA« 2 

-n x n||Aw 2 n 2 ± Av 2 - f(n) ) 

(36) 



with f(n) = ^v 2 A z 2 W(z). 

Note that for f = 0, Av 2 J. 0, 7 = (7) J. f , v% -> tt c /cj p 
and consequently W(z) — » z~ 2 , 136H reduces to the low- 
frequency limit of the cold plasma response tensor (lM2t with 
£ = 0. 

In the remainder of this paper, we will investigate 
whether the fast MSW eigenmode of \YA'2]i -induced by the 
GW- can scatter into an eigenmode of I36t when it propa- 
gates through a relativistic plasma region. 

5.3 Non-linear response tensor 

The quadratic response tensor for a relativistic magnetised 
pair plasma is obtained in the forward-scattering approach 
by averaging the second order ( in the field) cu rrent over an 
arbitrary distribution function l|Melrosell200lf) . We take the 
small gyro-radius limit to find: 



4m 2 



(37) 



x G*°"" (k , u)G* p " (ki , u)G* lp (k 2 , u) x 

ki k% 

T aa (k^u)T l 3 1 (k 2 \\u) + - — -T CT „(A;o||w)T ; 8 7 (fc 1 ||U)+ 



ko«u 

Kg 



ki\\u 
fcf 



k \\u 
k" 

p(kx\\u)T a -t{k2\\u) + T a f3(k2\\u)T aJ {k 1 \\u) + 



1 {k 2 \\u)Tp a {ko\\u) + 



k 2 \\u 
fcf 



fc 2 || U " ,v z " ' ' 1 k n u 

The cold plasma limit can be derived from 13711 by choosing 
F±(p) = n±5 4 (p — mu) (where again u = [1, 0] in the plasma 
rest frame) . The only components in 13 7H that depend on the 
sign of the charge are the gyrotropic terms which are pro- 
portional to r\ = q/\q\. In a pure pair plasma with equal dis- 
tributions for the electrons and positrons the non-gyrotropic 
terms are therefore the only ones that do not vanish in the 
sum over the species (see also the next Section EJ. 



6 SCATTERING 

As mentioned previously, the emission probability l|13t is 
very general and applies to any emission process due to an 



arbitrary source expressed as an external current density. 
The specific probability of scattering of a specific mode M, 
say a mhd wave, on an electron into a different mode P can 
now be calculated by inserting the extraneous single particle 
current 1261 1 into 1131 1 and interpreting the resulting emission 
as the scattered wave: 



WMp(k,ki,u) = 



3o- t (47t) 2 R M (k)R P (ki) 



7 



(38a 



(38b) 



'(fcjw , kip, ku — k\p)u p 



uj M {k)u}p(k\) 

x \a,Mp{k, k\,u)\ 2 5 [{k M \\ —ki P \\)t 

a,Mp(k,ki,u) = e* M ^(fyepvfa) { 

G al1 (k M ,u)r al 3(k M ^u)G* l3 ' / (kip , u) 

2m 
— c 

q 

where we have used [5{{k - k^u)] 2 = T5({k-k 1 )u)/(2n-y). 
We have included the non-linear scattering probability 
a M " p (fe.ivf, kip, ku — kip) due to the scattering of the mhd 
wave on the self-consistent field associated with the zeroth- 
order current related to the beat between the mhd wave 
and the scattered wave. As can be seen from l|38^ the total 
scattering probability is proportional to or oc q 4 /m 2 . The 
non-linear scattering probability has an additional (small) 
factor 2m/ q as can be seen from l|38bfl and consequently 
the gyrotropic terms are the only terms that do not van- 
ish when the contributions of the two species are added for 
a pure pair plasma. The components of the quadratic re- 
sponse, however, are proportional to (lj/Q c ) 4 and are negli- 
gible i n our long- wavelength approximation (see also lMel ros d 
l|l997MQedalin fc Machabelil ljl983h : lMikhailovskiil (jl98nfr T 

6.1 Polarization vectors 

The polarisation properties of the incident mhd wave and 
the scattered modes are determined by the eigenvectors of 
the wave operator A 11 " (k) as defined by 10 and Q with the 
linear response tensors derived in 132H and I36H . When the z- 
axis is chosen along the magnetic field and the wave vectors 
in the x — z plane, then the spatial parts of both the wave 
operator for the cold plasma and that for the relativistic 
an-isotropic plasma are of the form: 



An 





A13 





A 22 





A13 





A33 



(39) 



So the dispersion relations for the different wave modes, used 
in IIOL follow from: 

|Ay(fc)| = A 2 2(AnA33-A 2 3 ) = , (40) 

and the polarisation vectors are given by: 

et = {0,1,0}, (41a) 

={-Aia,0,Aii} , (41b) 



V(Ai 3 ) 2 + (An) 5 
1 

v / (Ai3) 2 + (A 11 ) : 



{An,0,Ai 3 } 



(41c) 



Clearly, e t always corresponds to a mode where the electric 
field is purely transverse to the magnetic field with disper- 
sion relation A22 = 0, such as the magneto-acoustic waves 
excited by a gravitational wave. 
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with 



Figure 4. Scattering geometry in observer (top) and comoving 
(bottom) frames. 



The modes e± that have an electric field component 
along the magnetic field correspond to the eigenfrequency 
solutions of 



(An - A 33 ) =F V , (2Ai3) 2 + (A 1 i-A 33 ) 2 = . (42) 



6.2 Response to an incident magneto-acoustic 
wave 

To evaluate the scattering probability Ij38^ explicitly, it is 
most convenient to oriented our coordinate axes such that 
the background magnetic field is aligned with the z-axis and 
the wave vector of the scattered wave (kp) lies in the x-z 
plane. The direction of the incident magnetosonic wave is 
then completely arbitrary, but we know from l|39p and l|41a|l 
that it is polarised perpendicular to both the magnetic field 
and the wave vector. Hence, we have 

km = {u, k sin 9 cos ip, k sin #sin ip, k cos 9} , (43a) 
e.M = {0, — sin ip, cos^, 0} , (43b) 
kp — {(jJi, k\ sin 0i, 0, fci cos 9i} , (43c) 

where 9 (9i) is angle of the incident (scattered) wave vector 
with respect to the magnetic field and ip (ipi) the angle with 
the a;- axis. 

Similarly, the polarisation of the scattered modes is 
found from the eigenvectors of the wave operator @ with 
the relativistic linear response tensor given by l|36[l using 



C2J: 

e t (fei) 
e_(fci) 

e+(fci) 



{0,1,0} , 
1 



VT+s 2 
1 



5f {-2,0,1} 



vTTs 



^{i 5 o,s} 



(44a) 
(44b) 

(44c) 



„ k\u\ sin 9\ cos t 
cu'f — k\u\ cos 2 



(45) 



With this choice of geometry, the explicit expression for 
the linear scattering response summed over the two species 
such that the gyrotropic terms vanish, is: 

a M "(fe, ki,u) = G° !fl (A;M,it)r Q( 3(fc M ||M)G* /3l '(feip,M) , 



with the spatial components written explicitly as: 

(ku) 2 



1 



sin 6 cos ifj 
ku 



aij(k,ki,u) 


1 

7/3fc sin 8 sin ip 
ku 



(ku) 2 - Q 2 



"Y0ki sin 6\ 
ku 





(46) 



and * = 4^^ < fc ;> ~?r + (1 - 7 2 ) 



2 \ fc sin cos i/> ki sin 0i 
ku k\u 



6.3 Probability of transverse-transverse scattering 

The probability that the incident (transverse) magneto- 
acoustic wave scatters into a transversely polarised electro- 
magnetic wave, as in (I44at . is proportional to (the square 
of): 



aMt{k,ki,u) = eA/ fl (fc)e tI/ (fci)a' ly (fc,A;i,'u) 
(ku) 2 _ 1 



(47) 



(ku) 



-n 2 



7 2 (^-fc / 3cose) 2 



In the two limits of nearly parallel and perpendicular prop- 
agation and assuming that 70; <C fi c , this is approximately 
~ (Luf-yQc) 2 and ~ (-yuj/Q c ) 2 , respectively. 

Furthermore, it can be easily seen from the linear re- 
sponse tensors in i'A'21 and l|36^ , with the polarisation vectors 
of the corresponding transverse eigenmodes given in I|43b0 
and 144a> . respectively, that the amplitudes 



a.M,t(k) = e k M>ta (k)e M ,tv(k)ct til/ (k M ,t) \ u 



(48) 



are independent of u) for the transverse magneto-acoustic 
and electromagnetic wave modes. Hence, the ratio of electric 
to total energy given by p6a^ is Rm = Rt = 5 just like for 
transverse modes in an isotropic non-dispersive plasma. 

Compared to the probability of scattering of a fast 
magneto-acoustic mode to a transverse electromagnetic 
mode in unmagnetised vacuum, this probability is sup- 
pressed in a strongly magnetised plasma by the factor in 

squared. 

6.4 Probability of transverse-LO scattering 

The probability that the transverse magneto-acoustic mode 
is scattered into a quasi-transverse Langmuir-ordinary mode 
is obtained by contracting l|46|l with the polarisation vectors 
l|43bfl and either l|44bfl or l|44cfl . 

The only component of 146H . however, that is of lower 
order in 1/Q. C is the purely longitudinal zz component, which 
has a term independent of tt c . This component drops out, 
though, as soon as (I46H is contracted with the polarisation 
vector of the purely transverse incident magnetosonic wave. 
We find that the probability for scattering into a LO mode 
is equal or less than 1471 depending on the relevant angles, 
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the maximum being when the incident wave is perpendicular 
to the magnetic wave and the electric field of the scattered 
wave lies in the plane spanned by the electric field of the 
magneto-acoustic wave and the magnetic field. 



6.5 Emitted power 

In the previous section we derived the probability I38H that 
one quantum in the MSW mode scatters into one quantum 
in another mode. To calculate the actual rate at which scat- 
tered waves are emitted we have to average 13811 over the 
distribution of particles. One subtlety is that the scatter- 
ing probability per unit time is clearly a frame-dependent 
quantity. We should therefore integrate either the scattering 
probability r jWMp{k,ki,u) per unit proper time along the 
particles motion, multiplied by the proper number density 
d 4 pF(p) or equivalently the probability per unit time over 
the physical number density 'yd 4 pF(p), similar to the aver- 
aging of the response tensor over the distribution function 
performed in 1281 . 

Appealing to the semi-classical formalism of detailed 
balance, we find a kinetic equation for the increase of wave 
energy in the scattered modes: 



DWpjk! 
Dt 



(2^)3 



dp\\WMp{k, ki,u)g{jp\\ 



x I —W M (k) - W P (ki 



(49) 



where we only consider spontaneous scattering of M — ► P 
and P —y M and neglect induced scattering 6 . 

The total power in 1491 is calculated for the scattering 
on the relativistic distribution of particles along the mag- 
netic field in the frame that is comoving with the average 
bulk flow carrying the magnetic field. Because the power is 
Lorentz invariant (when the scattered radiation has front- 
back symmetry in the comoving frame) it is the same in the 
observer frame. 

However, because of propagation and angular effects 
there is a difference between the emitte d power P e and the re- 
ceive d power P r in the observer frame iRvbicki fc Lightmanl 
Il979l) . The scattered radiation is emitted during the time it 
takes the GW and MSW to cross a scale-height, i.e. Ti/c, but 
because this shell propagates at a relativistic velocity f3, the 
radiation is received in an interval Ti/^-ypc) with a power 



(50) 



These relativistic effects also play an important role 
in the physics o f gam ma-ray bursts as reviewed in 
IZhang &: Meszarosl (2004). Hence the time-integrated emit- 
ted and received energies are the same but (with sufficient 
time resolution) the observed flux is much higher. 



7 APPLICATION TO NS-NS MERGER 

We consider the specific environment of a highly rela- 
tivistic electron-positron wind surrounding the merger of 



6 D/Dt is a time-derivative operator that allows for a possible 
slow variation of the m edium itself as D/Dt = (dwM /dk^dp + 
{du M / 9x^)0 /dki 1 as in lMelrosd bOOll) . 



two neutron stars as illustrated in Fig. [31 and apply the 
general results for the scattering of low-frequency fast 
msws into escaping electromagnetic waves as derived in 
the pr e vious sections f r om th e framework inlGedalin et alJ 
Jl998ft: IMelrose et alJ dl999l): IMelrose fc Gedalinl <199Sft: 
Gcdal in et al .1 l|200lft : lLuo fc Melrose! <|200lh : lOedalin et alJ 

We assume that the in the last phase of coalescence, 
the orbital motion dominates the spin of the individual neu- 
tron stars. The orbital angular frequency Q.^ then deter- 
mines the radius of the light-cylinder R\ c = c/fif, where the 
magnetic field lines that are anchored on the stellar polar 
caps break open because co-rotation of the plasma which 
is frozen onto the field lines would result in an unlimited 
increase in inertia if the filed lines were closed at this dis- 
tance. For a more extensive treatment of how the spinning 
magnetic dipole of a neutron star sets up a voltage jump 
and an electric field at the stellar polar caps and extracts 
charges fro m the neutron star fillin g the magnetosphere, see 
for instance iGurevich et al] <ll993T) . Suffice to say that dur- 
ing the last stage of spiral-in, the double neutron star is 
assumed to behave essentially as a single millisecond pul- 
sar with a charge separated closed magnetosphere and an 
ultra-relativistic force-free electron-positron plasma wind. 

An important observation is that the relativistic lep- 
tonic wind is already present before the actual merger. It 
doesn't have to be generated by the merger event, which 
is a problematic issue in collapsar models for gamma-ray 
bursts, where inside a collapsing massive star a region has 
to be evacuated of baryons and an almost mass-less jet has to 
form and then pierce through the entire star without slowing 
down much or losing its collimation. 

As the neutron stars finally merge, a large fraction of the 
binding energy of the entire system is released in the form of 
gravitational waves that propagate through the pre-existing 
strongly magnetised plasma wind. Since a magnetised per- 
fect fluid possesses an an-isotropic pressure due to the mag- 
netic stress, the GW couples to the plasma and excites mag- 
netohydrodynamic waves. In Paper II we found that Gws 
can excite both Alfven and (slow and fast) magneto-acoustic 
waves but that coupling to the fast mode is the most efficient 
because the dispersion relations are nearly the same in a ten- 
uous, relativistic magneto-plasma and coherent interaction 
is always possible. 

Even though we proved that the energy transfer from 
GW to MSW in the plasma can be substantial, we still have 
to find an efficient radiation process to release this energy. 
The MSW cannot escape the plasma directly, so we will now 
investigate whether the MSW are likely to inverse Compton 
scatter into electromagnetic radiation in the radio regime. 

7.1 MSW amplitude and energy 

The excitation of fast magneto-acoustic waves by a GW prop- 
agating through a uniform magnetic field was discussed in 
Section [5] The orientation of the ambient magnetic field is 
now chosen in the z-direction such that the wind flows out in 
the ^-direction with average bulk velocity f3 and correspond- 
ing Lorentz factor yp. The wave vector of the gravitational 
waves, k gw , is allowed to have an arbitrary angle 8 with 
respect to the magnetic field (and the wave vector of the 
excited MSW will have the same orientation). In reality, the 
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GW are emitted more or less radially and k gw will be almost 
perpendicular to the toroidal magnetic field in the wind in 
the observer frame (i.e. as in the top Fig. |1J. 

In the approximation of a uniform background magnetic 
field, the amplitude and energy of the excited MSW are given 
in the observer frame by 

B\ ' (x) ~ ■ — ^ — -smo k gvi x e B y ' , (51a) 

27/3 

w w = [ d 3 x l^f . ( 5i b) 

Jv 2 Mo 

W M = T 00 ~ W (E) + W (B) ~ 2W [B) . (51c) 

This result will be generalised to a non-uniform background 
in the following sections. 




Figure 5. Wind collimation by the magnetic field for k = 1 (left) 
and k = i (right). 



7.2 Eikonal approximation 

We generalised (I51at to a non-uniform magnetoplasma wind 
by appealing to the eikonal or geometric optics two-length- 
scale expansion (also known as the wkb approximation). 
When the wavelength of the modes under consideration 
is small compared to the scale of spatial variation in the 
background 7 , the waves can still be considered locally pla- 
nar and monochromatic. Furthermore, in the Poynting flux 
dominated regime where the gas pressure is negligible (Sec- 
tion the classical Alfven velocity va oc B/i/n is con- 
stant since B ~ £?0 oc 1/7? well outside the light-cylinder 
and n oc B z oc 1/R 2 (which will be derived in ( I53H ). 
Consequently, the relativistic Alfven phase velocity 1/ma = 
l + l/v A is also constant and the dispersion relation cu ~ kuA 
is spatially non-dispersive even in the inhomogenous back- 
ground. 

The important space-dependent background quantities 
are the ambient magnetic field S^(r) and plasma density 
n(r), and the corresponding gyro-frequency f2 c (r) oc S'°'(r) 
and plasma frequency u) p i(r) oc \J n(r). Note that the grav- 
itational waves h+ (r) are also treated in the 
geometric optics approximation as a field living on a flat 
space-time background, neglecting the curvature due to the 
neutron star and the (largely magnetic) energy density in 
the wind. The GW amplitude is estimated from the total en- 
ergy released by the merger in GW at a distance from the 
source where we are in the far field approximation, such that 
the GW can be treated as linear, plane wave space-time per- 
turbations. The GW amplitude then falls off as 1/r and we 
can scale it such that h+(r) — hi c R\ c /r. 

Finally, the linear growth of the excited MSW pertur- 
bation in l)51a|l can be generalised to the inhomogenous jet 
by considering growth over single scale heights TL (see I55H ). 
Putting all of this together, we find 

B< 4> } ( r ) ^ hi c k gw Ri c c ik x „(r-t) ( 52a) 

= ^ /e lfcsw(r - t) , (52b) 



where the interaction factor ^/ is now a constant deter- 
mined by 152aH and the dependence of B^\ B^ on r will 
be determined in the next sections. 

Note that because the GW and the MSW both have es- 
sentially the same dispersion relation, and the energy flux, 
S'ic_B^ lhd /(47r), through each surface Si (as in Fig. |5J is 
a constant because of l|52b|l . both waves can propagate in 
phase and interact coherently over a long distance even in 
the non-uniform background. 

7.3 Magnetic field configuration 

The morphology of the 'pulsar' wind in which the GW-MSW 
interaction and the subsequent inverse Compton scattering 
take place, is fully determined by the magnetic field geom- 
etry since the particles are 'frozen' to the field lines. We 
assume that the wind is steady and has a constant bulk 
velocity v x , the magnetic field is axially symmetric and 
force-free, and we use the fact that the field is solenoidal. 
For such a field, the toroidal component follows from the 
torque free condition a s RB^, = const, along each field line 
iLiist fc Schliiterlll954l) . where (R,<f>,x) denote cylindrical 
coordinates. Note that because of <- r >2t . the discussion in 
this section is valid for both and oc B° and we will 
omit the superscripts to avoid clutter. 

The poloidal magnetic field B p through any flux tube 
cross-section A5* must be proportional to the particle flux 
through AS because of flux conservation and the constant 
flow velocity v: B P AS = const., nvAS = const., and thus 
B p oc n. In particular, we will assume that the particle and 
(poloidal) magnetic flux through an azimuthal cross-section 
of the entire jet will approximately be conserved (as illus- 
trated in the right Fig. 0: B x Si = B x (irRf) ~ const. 

Summarising, we have 

£^0R)oci, ft,(fl,(a;))(xn(/%(x))(x , (53) 

where Rj(x) is the radius of the jet which is related to the 
opening angle by o (x) — 2 arctan Rj(x)/x. 



7 We assume a time-independent background, but one can ob- 
viously also use the eikonal approximation on a time-dependent 
background when the period of the waves is small compared to 
the time scale of variations in the background. 



7.4 Parametrised wind collimation 

In Papers I— III we have discussed different degrees of col- 
limation of the wind. We can generalise these results by 
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parametrising the collimation as in Fig. such that the ra- 
dius Rj of the jet as a function of the vertical distance x is 
given by Rj = x K . From l|530 then follow the more useful 
expressions for the vertical dependence 



Bd, oc x 



V = 



B x oc x 



tvR ax x x 



(54a) 
(54b) 



We can define a scale-height TC (as in Section 17,21 by the 
distance over which the toroidal magnetic field falls off to 
half its strength: 



B<t,{x + H) _ i 
B*(x) ~ 2 



H 



l)x 



(55) 



we find that for k = 1/2 the volume integrated power is a 
constant and that it increases with distance only for 2re > 1. 
Hence, we will next consider a collimation of re — 1/3. 



7.4.3 



Clearly, for re — 1/3 the magnetic field and hence the MSW 
energy density remain much higher throughout the wind, 
but the decrease in integration volume is also significant as 
can be seen from Ij54bfl and l|59JI. The total energy deposited 
in MSW by the GW for the same parameters as in Section lV.4.11 
and over the last scale-height, which is now Tt — 7r from 1551 
and with r — \/x 2 + R 2 ~ r is 



7.4.1 k = 1 

In Paper I we assum ed a wide conical wind with k = 1 (see 
also iKuiipersI l)200lfl 1 as depicted in Fig- E] (left) which has 



Bd> oc 



Wm oc i 



B x oc 



V oc r 



(56a) 
(56b) 



with the radial distance r = \J R 2 + x 2 ~ x. The vol- 
ume integrated energy increases linearly with the interaction 
length-scale set by the scale-height Ti = r from I55H . We will 
only consider the energy deposited in the last scale-height of 
the force-free wind, as the radiation will escape most easily 
there. 

For a magnetar strength surface magnetic field of ~ 
10 12 T, a millisecond orbit, a Lorentz factor of 7/3 ~ 100, a 
GW amplitude h\ c = 10 _3 7?i c /r and i? max ~ 10 14 m we find 
that the total energy transferred from the GW to the MSW 
in the plasma is 



Wtot(«=l) = 



7.4.2 K=\ 



W M (r)ivr 2 dr ~ 10 35 J . (57) 



In Paper III we made a numerical estimate of the power 
density in inverse Compton scattered radiation by neglect- 
ing the magnetic field in t he scattering process and using 
iRvbicki fc Lightmanl (Il979l) and 



Pic(r 



: C7|/3 2 a T n(r)W M (r) 



(58) 



with 7/3 and n the Lorentz factor and particle density of 
the scattering particles, respectively, and as before, Wm the 
energy density in incident MSW. 

The power in scattered radiation not only depends on 
the energy density in the mhd waves Wm(t) but also on the 
density of scattering particles n(r), and the parametrised 
distance dependence of the IC power density follows from 
lEU as 



-Pic oc x 



(59) 



To obtain a volume integrated power we assumed in Pa- 
per III a collimation parameter of re = 1/2 in which case 
the spatial dependence in 159H is the same 1/r 2 decay as for 
the MSW energy density for the re = 1 collimation in l|56bfl . 
However, we overestimated the total scattered power by in- 
tegrating over a spherical shell oc r 3 . From 154bt and 15911 



Wtrt(K=£) 



*/8 



WM(r)nr 3 dr < 



10 28 J 



(60) 



Note that from 1541 1 and 15211 it is clear that the volume 
integrated energy in the MSW always increases linearly with 
distance, which is a consequence of flux conservation and 
the discussion leading to 15211 (such t hat hTC = const.). 

If we repeat the estimate in iMoorteat fc Kuiipersl 

(2005) and integrate l|58|l but now for re = 1/3 such that 



n(r) 



4Me QbBi ( 



Ric 

r 



(61) 



derived from the Goldreich- Julian density in terms of the 
multiplicity M ~ 10 3 and orbital frequency Qb, and 



-Pic, tot = C7/3 0T 



*/8 



WAi(r)n(r)-Kr 3 dr > 



10 36 w 



(62) 

which for a binary merger close enough for a LIGO detection, 
e.g. in the Virgo cluster at a distance of ~ 1 Gpc, would yield 
a flux on earth observable by lofar. The duration of this 
kind of transient signal is relatively short: At = H/(2'fgC) w 
15 sec, but still much longer than the cosmic ray events that 
LOFAR will also detect. 

Unfortunately, we have learned in this paper that a 
more careful analysis including the polarisation of the in- 
cident and scattered wave modes with respect to the strong 
ambient magnetic field leads to a strong suppression of the 
emitted power in the scattered radio waves. We can find 
a relation similar to l|62|l from 1)49^ by considering mono- 
energetic particles with, say, a Lorentz factor 7/3 (/3) ~ 10 2 
and distribution function g(p\\) = <5(pii — mf3)n/(2nm). To 
evaluate the integral over wavelengths in l|49|l we assume 
that the incident GW and mhd waves are mono-chromatic 
with k = feg W . In comparing I62H with 149H . however, the 
latter has an additional factor w, the scattering probability, 
which depends on the small value of 1471 1. 

In Fig. [(J we plot the (non-relativistic) gyro-frequency 
Q c for a pulsar strength magnetic field (B* = 10 8 T) as 
a function of the distance in the wind, together with the 
plasma frequency (middle line) and for comparison a char- 
acteristic LOFAR frequency of 100 MHz (bottom line). Al- 
though the relativistic gyro- and plasma- frequencies are sup- 
pressed by the Lorentz factor, it is clear from Fig. |S| that 
even for 7 ~ 10 2 -10 7 , the gyro-frequency greatly exceeds 
the MSW frequency and even the frequency of the scattered 
waves everywhere in the wind. 
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Figure 6. Relevant frequencies as a function of distance. 
Top: (non-relativistic) gyro-frequency; Middle: (non-relativistic) 
plasma frequency; Bottom: characteristic LOFAR frequency of 
100 MHz. 



7.5 Synchrotron radiation 

In the discussion so far we have assumed that the ambi- 
ent magnetic field in the wind is sufficiently strong that 
any transverse momentum is radiated on a time-scale short 
to that of the other processes considered. This allowed us 
to treat the particle distribution function as strictly one- 
dimensional. However, if the magnetic energy density drops 
sufficiently, the small gyro-radius approximation in Sec- 
tion l4.2l is not valid anymore. Particles acquire a pitch-angle 
and synchrotron radiation can become important. 

We calculate whether, as the magnetic field drops to 
its minimum strength, either the radius of curvature of the 
magnetic field or the MSW wavelength could become com- 
parable to the maximum gyro-radius r c (i? max ) of particles 
that acquire some transverse momentum. In the former case, 
the adiabatic invariant of the particle motion, the magnetic 
flux through the particles orbit nr^B — const, could be 
violated. 

However, we find that even when the particles orbit the 
magnetic field with relativistic velocity v± ~ Vu ~ c and 
corresponding Lorentz factor 7 C ~ 100, we find 



£^0.3 lO-VM 

n c LiooJ 



B, 



10 s T 



(63) 



which is clearly much smaller than the curvature radius of 
the magnetic field (which is proportional to TL 1 ^) and the 
MSW wavelength A msw = 2nc/ui sv , — 150 km. 

The smallness of r c is partly due to the rapid rota- 
tion of the source, causing the magnetic dipole field to 
drop only by a factor ~ 125 from the surface to the light- 
cylinder. In the wind, the field remains strong because of 
the force-free condition combined with the collimation. For 
k = 1/3, it only decreases by another factor ~ 10 3 in 
the entire wind and at -R max = 10 14 m the field is still 
0.7 10 3 T. Even for k = 1, though, we still 



Bf\Rmzx) ' 

have B ( ° } (R 



0.4 10" 3 T and r c {R n 



0.4 km 



(and now the radius of curvature of the magnetic field is 
ccHcc -R max ). 

The comparison of r c to the wavelength of the magne- 
tosonic waves is related to the factor 147H suppressing the 



inverse-Compton scattering 
r c {r) _ 



2tt fi c (r) ' ^ 

Consequently, if the electrons and positrons are initially con- 
fined to the magnetic field lines, it is unlikely that they will 
acquire transverse momentum and emit synchrotron radia- 
tion anywhere in the wind. 



8 DISCUSSION 

Since a strong magnetic field is essential to have efficient 
interaction between GW and MSW, and millisecond rotation 
(with a correspondingly tight light- cylinder at R\ c ~ 50 km) 
is necessary to have sufficiently large frequencies for the MSW 
and the scattered waves, 'yui/Q will be small throughout the 
wind. For the range of parameters studied here, frequencies 
in the LOFAR band never exceed the local relativistic gyro- 
frequency within the entire extent of the wind. 

Furthermore, we know from Paper I and Paper II that 
the GW most efficiently excite purely transverse magneto- 
acoustic waves, so the scattering channel from a (quasi- 
longitudinal to another (quasi-)longitudinal mode, which 
doesn't suffer from cyclotron suppression, is not a viable al- 
ternative. The reason that scattering is suppressed is that 
the incident MSW does not have an electric field compo- 
nent along the magnetic field, and since the electrons and 
positrons are assumed to be constrained along the field lines 
they are unable to oscillate in the applied electric field per- 
turbation. 

Given the large amount of energy transferred from the 
GW to the MSW in the plasma (1571 it is very likely that some 
radiation mechanism will produce observable radiation in a 
higher frequency band. The analysis presented here proves 
that the simplest possibilities, synchrotron radiation, and 
inverse Compton scattering of the GW induced MSW into 
low- frequency radio waves are ineffective. We have yet to 
find a more viable alternative. 



9 CONCLUSIONS 

In this paper we have studied scattering processes in an 
essentially one-dimensional, intrinsically relativistic, mag- 
netised plasma in some detail with the aim of investigat- 
ing whether very low-frequency magnetosonic waves ex- 
cited by a passing GW can scatter on the relativistic elec- 
trons and positrons in the wind to produce radio emis- 
sion in the frequency regime observable with LOFAR. To 
treat the highly relativistic system properly, we used the 
covariant and gauge independent t heory of plasma dy- 
namic s developed in lOedalin et all (ll99St): IMelrose et alJ 



il999ft: IMelrose fe Gedalinl (I1999T): iGedalin et alJ J200_. . 
iLuo fc Melrose! l|200lft : lGedalin et alJ l)2002l) . This approach 
also allows clear comparison with our general relativistic 
mhd treatment in Paper I and Paper II of the propagation of 
gravitational waves through a magnetised plasma. By spec- 
ifying to the temporal gauge, our results are translated to 
more familiar vector expressions. To our knowledge, the ex- 
plicit expressions for the scattering probabilities for inverse 
Compton scattering in a strongly magnetised, force-free and 
intrinsically relativistic plasma jet have not been derived 
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before in the literature. These results are useful as well in 
the theory of ordinary pulsar winds and the jets of X-ray 
binaries and active galactic nuclei. 

We assumed that the magneto-acoustic waves are ex- 
cited in a cold, essentially mono-energetic plasma that sat- 
isfies the ideal mhd condition. As they propagate outwards 
these waves are assumed to encounter either high energy 
tails of the particle distribution or, in general, a region with 
an intrinsically relativistic particle distribution on which 
they scatter. Since the plasma is strongly magnetised, all 
momentum transverse to the magnetic field is synchrotron 
radiated on a short time-scale, and as a result the charged 
particles are strictly confined to the magnetic field lines. In 
evaluating the response of the plasma, we therefore average 
over the relativistic one-dimensional Jiittner-Synge distribu- 
tion function. 

To simplify our results we restrict ourselves to the long 
wavelength regime with ui <C £l c , which is justified in Sec- 
tion[7| Since the incident magnetosonic waves excited by the 
GW are strictly transverse to the magnetic field, the scatter- 
ing probability for all possible scattered modes is found to 
be suppressed by a factor (jLu/fl c ) 4 . Thus, inverse Comp- 
ton scattering of the kHz MSW to ~ 100 MHz radio waves 
does not seem to be a preferred radiation mechanism for the 
energy dumped in the plasma by the GW. 

The covariant formalism for relativistic plasma dynam- 
ics that we have used was developed primarily to study emis- 
sion mechanisms in pulsar magnetospheres where the dipo- 
lar field lines point out almost radially from the polar caps, 
the charged particles flow out along these field lines and the 
most interesting scattering processes usually have both the 
incident and scattered wave vectors almost parallel to the 
magnetic field. The situation that we have studied is the op- 
posite: in the force-free wind well outside the light cylinder, 
the magnetic field is dominated by the toroidal component 
and is essentially perpendicular to the incident and scattered 
modes of interest (in the observer frame). For a millisecond 
pulsar and a slightly collimated wind, the magnetic field 
is very strong throughout the wind and the particles move 
along the field lines. Note, however, that the particles have 
relativistic velocities both in the radial direction and along 
the magnetic field lines and in the laboratory frame appear 
to be moving outwards like 'beads on a spiralling wire' pre- 
dominantly in the radial direction. Consequently, the angle 
between both the incident and scattered wave vectors and 
the velocity of the particles along the magnetic field can be 
quite small in the lab frame (due to relativistic beaming), 
while they are almost at right angles in the (average) rest 
frame. 
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